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ABSTRACT 

First-generation (Population III) stars in the universe play an important role 
in early enrichment of heavy elements in galaxies and intergalactic medium and 
thus affect the history of galaxies. The physical and chemical properties of pri- 
mordial gas clouds are significantly different from those of present-day gas clouds 
observed in the nearby universe because the primordial gas clouds do not contain 
any heavy elements which are important coolants in the gas. Previous theoretical 
considerations have suggested that typical masses of the first-generation stars are 
between several M & and « 1OM although it has been argued that the formation 
of very massive stars (e.g., > 1OOM ) is also likely. If stars with several M are 
most popular ones at the epoch of galaxy formation, most stars will evolve to 
hot (e.g., > 10 5 K), luminous (~ 1O 4 L ) stars with gaseous and dusty envelope 
prior to going to die as white dwarf stars. Although the duration of this phase is 
short (e.g., ~ 10 5 yr), such evolved stars could contribute both to the ionization 
of gas in galaxies and to the production of a lot of dust grains if the formation of 
intermediate-mass stars is highly enhanced. We compare gaseous emission-line 
properties of such nebulae with some interesting high-redshift galaxies such as 
IRAS F10214+4724 and powerful radio galaxies. 

Subject headings: galaxies: formation - galaxies: evolution - galaxies: nuclei - 
galaxies: individual (IRAS F10214+4724) - galaxies: radio sources - planetary 
nebulae 
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1. INTRODUCTION 

In the standard big bang cosmology, the amount of heavy elements synthesized in the 
universe until the decoupling of matter from primordial radiation is negligibly small; i.e., 
Z < 10~ 10 where Z is the mass fraction of heavy elements (Wagoner, Fowler, & Hoyle 1967). 
On the other hand, gas clouds in the present-day galaxies have an appreciable amount of 
heavy elements Z ~ 0.01; e.g., the solar heavy element abundance is Z Q « 0.02. This means 
that chemical and thermal properties of gas clouds in primeval galaxies at high redshift are 
significantly different from those in the present-day galaxies. 

Just after the recombination of the universe, the hydrogen gas could be in atomic phase 
at a temperature of 10 4 K. Then the formation of primordial hydrogen molecules can 
proceed through the gas phase reaction (e.g., Peebles & Dicke 1968). Since the rotational 
and vibrational transitions in hydrogen molecules can cool the gas further, the first generation 
stars can be born in such gas clouds at redshift ~ 10 - 30. Since the temperature is still much 
higher than a typical kinetic temperature of molecular gas clouds in the present-day galaxies 
(e.g., ~ 10 K), the initial mass function (IMF) in primeval galaxies may be different from that 
in the present-day galaxies; i.e., a Salpeter-like IMF (Salpeter 1955; Scalo 1986). Therefore, 
one of fundamental questions on the star formation in primeval galaxies is; "What kind of 
stars (i.e., Population III stars) were made in such extremely metal deficient gas clouds ?" 
[Matsuda, Sato, & Takeda 1969; Yoneyama 1972; Yoshii & Saio 1986 (hereafter YS86); Silk 
1977a; Tegmark et al. 1997; Nakamura & Umemura 1999, 2001a, 2001b; Bromm et al. 1999, 
2001, 2002; Coppi et al. 2001]. 

Applying the opacity-limited fragmentation theory [the fragmentation of gas clouds due 
to the gravitational instability continues until the fragments become to be optically thick 
to the cooling radiation (Silk 1977b, 1977c)], YS86 first derived an IMF for a hydrogen- 
helium gas cloud. They found that typical masses of the first-generation stars are from 
several M & to 10 M Q . Nakamura & Umemura (1999) obtained similar results based on 
their one-dimensional hydrodynamical calculations together with nonequilibrium processes 
for hydrogen molecule formation. Abia et al. (2001) pointed out that the very large C and N 
enhancements observed in the extremely metal-poor stars in our Galaxy favor a Population 
III IMF peaked at intermediate-mass stars. 

Recently, on one hand, some groups investigated the collapse and fragmentation of 
primordial, metal-free gas by numerical simulations (Abel, Bryan, & Norman 2000, 2002; 
Bromm et al. 1999, 2001, 2002; Coppi et al. 2001). They showed that the first generation 
gaseous clumps have masses of > 1000M o which corresponds to the Jeans mass of the gas 
with temperature T ~ a few 100 K and the density n ~ 10 3 — 10 4 cm~ 3 . As Bromm et 
al. (2002) assumed, if the clump masses are indicative of the final stellar mass, Population 
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III stars may be more massive than 1OOM . Nakamura & Umemura (2001a, 2001b) have 
developed their hydrodynamical simulations and found that the initial mass function of 
Population III stars may be bimodal with peaks of ~ 1 - 1OM and ps 1OOM . The 
mass depends on the initial density and the initial fraction of molecular hydrogen (xh 2 ,o) 
(Nakamura & Umemura 2001b). If :th 2 ,o is larger than ~ 3 x 10~ 3 for low-density (n < 10 5 
cm -3 ), the effective HD cooling makes the high mass peak small to ~ 1OM . However, if 
Xu 2 ,o is smaller than ~ 3 x 10~ 3 , very massive stars with ~ 100M o are formed preferentially. 

In this way, the mass of first-generation stars is controversial. It seems also hard to 
examine observationally which mode is more feasible; i.e., Population III stars are either 
intermediate-mass star or very massive ones. If most Population III stars are as massive 
as ~ 100M o , they could explode quickly and then contribute to the chemical enrichment 
of either galaxies or intergalactic medium or both (e.g., Madau, Ferrara, & Rees 2001 and 
references therein). Since usual Population II stars could be formed in cold gas clouds 
polluted chemically by first Population III stars, we cannot examine observationally whether 
or not such very massive stars were really formed in the primordial gas clouds. Furthermore, 
according to Nakamura & Umemura (2001a), proto-galactic clouds with a higher density tend 
to form intermediate- mass stars. If such proto-galactic systems are present, their evolution 
can be governed by evolution of the intermediate-mass stars whose lifetime is ~ 10 8 yr, 
being much longer than those of very massive stars. We investigate what happens and 
follows in such bursts of formation of intermediate-mass stars and discuss effects on the 
further evolution of galaxies, adopting the IMF of YS86. In particular, intermediate- mass 
stars experience inevitably the planetary-nebula-nucleus (PNN) phase, providing a lot of 
high-energy photons (e.g., Vassiliadis & Wood 1994). Therefore, it is expected that the 
galaxies appear to be luminous emission-line galaxies during the PNN phase. We discuss 
how they look and compare our model results with observations of some interesting high- 
redshift galaxies. 

2. EVOLUTION OF POPULATION-III STARBURSTS IN FORMING 

GALAXIES 

2.1. Initial Mass Function for Population III Stars derived by YS86 

In many previous galaxy evolution models, it has been often assumed that stars are 
formed with a Salpeter-like initial mass function (IMF) in a galaxy; 



4>(m) = f3m ■ 



(1) 
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where m is the stellar mass in units of M and (3 is a normalization constant determined by 
the relation 



'mi 

which leads to 



/ <f)(m)dm = 1, (2) 



(/j - ik- 1 m 

In this formulation, there are three free parameters; the power index (//), and the upper 
and lower mass limits of the IMF (m u and mi). In order to describe the evolution of the 
solar neighborhood, one adopts conventionally ji = 1.35, mi = 0.1M Q , and m u = 60M Q or 
1OOM . However, it is highly uncertain whether or not this Salpeter-like IMF is applicable 
to the star formation history in primeval galaxies. 

YS86 investigated the IMF for a hydrogen-helium gas cloud in detail based on the 
opacity-limited fragmentation theory. They found that typical masses of the first-generation 
stars are from several M & to 10 M & , depending on the assumed mass-luminosity relation for 
protostars derived by Silk (1977c); assuming that the opacity sources within the protostar 
is due to grains, he obtained the protostar luminosity 



L(m) = el0 2A7 (m/M Q r A=>, (4) 

where e = (Z/Zq)^ 1 and rj = 3. Since this equation cannot be applicable to the luminosity 
of stars with Z < 10~ 4 Z & , YS86 treated e as well as i] as free parameters in the ranges 
!0 < e < 100 and 1.5 < 77 < 3. Then YS86 found that the mass at the peak of the 
IMF is ~ 4M Q for i] = 3 and 10M q for i] = 1.5. Another important conclusion in 
YS86 is that stars more massive than these peaks are formed less efficiently than those 
expected from the Salpeter IMF with a = 1.35. In addition, stars less massive than these 
peaks are also formed less efficiently than those expected from the same Salpeter IMF (see 
Figure 3 in YS86). Therefore, the IMFs derived by YS86 are significantly different from 
the Salpeter IMFs. We therefore adopt the IMFs of YS86 and investigate luminosity and 
spectral evolution of galaxies for various sets of the parameters. 

For convenience, we present a parametric form of the YS86 IMF in the mass range 
between 1 M and 60 M for r] = 3 (YS86a), 
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4>{m) = < 



' 1.5944 x 1(T 2 m 2 - 
4.1728 x 1(T 2 m - 95 
1.4462 x KT 1 m° 
1.0813 m- 1 - 25 



I < m/M < 2.5 
2.5 < m/M & < 3.7 
3.7 < m/M & < 5 

5 < m/M Q < 7.5 
7.5 < m/M < 11 

II < m/M Q < 60 



(5) 



2.0081 x 10 mr 2 - 1 
2.2088 x 10 2 m~ 3 - 7 



and that for 77 = 1.5 (YS86b), 



' 7.3036 x 10~ 4 m 3 - 5 
4.4629 x 10~ 3 m 1 - 8 



1 < m/M Q < 2.9 
2.9 < m/M Q < 5 
5 < m/M Q < 8 
8 < m/M Q < 12 
12 < m/M G < 60 



4>{m) = < 8.0865 x 10" 2 m° 
1.2849 m- 1 - 33 



(6) 



1.4671 x 10 m 



-2.31 



Here the above numerical coefficients are determined by the following normalization for each 
case, 



In Figure 1, we show the shapes of IMF which we use in this paper. Since the massive end 
slope of YS86a IMF is steepest, the number of massive stars (> 2OM ) is smallest among 
the three, although the number of intermediate-mass star exceeds that of Salpeter IMF. 



Prior going to investigating the luminosity and spectral evolution of galaxies using the 
YS86 IMFs given in equations (5) and (6), we estimate how populous intermediate- mass 
stars are in the YS86 IMFs with respect to those expected from usual Salpeter-like IMFs. 
Since it is expected that the most spectacular effect of intermediate-mass enhanced star 
formation is brought by a number of luminous PNNs with very high temperature ~ 10 6 K 
(e.g., Vassiliadis & Wood 1994). Since the mass of such high-temperature PNN progenitors 
lies in a range between 4M and 5M , we estimate how many stars within this mass range 
are formed. 

If we adopt a Salpeter IMF, the number of stars with a mass range m\ < < m 2 
formed in the gas with 1 M is estimated as 




^60 



(7) 



2.2. How Populous are Intermediate-Mass Stars ? 
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Taking m x = 4M Q and m 2 = 5M , we obtain the numbers for the following two cases; 1) /i 
= 1.35, mi = O.1M , and m u = 6OM , and 2) // = 1.35, mi = 1M , and m u = 6OM (see 
Table 1). 

If we adopt a YS86 IMF, we can estimate the number of stars with a mass range 
4M < < 5M by integrating equations (5) and (6) for 77 = 1.5 and rj — 3, respectively. 
The results are shown in Table 1. We also show the number if all stars have a mass of 5 M ; 
i.e., the IMF is described as a 5 function with m* = 5M . 

As shown in Table 1, the number of intermediate-mass stars with m* = 4M - 5M for 
the YS86 IMF with rj = 1.5 is nearly the same as that for the Salpeter IMF with mi = 1M . 
However, that for the YS86 IMF with 77 = 3 is larger by a factor of 2.4 than that for the 
Salpeter IMF with mi = 1M . For both cases, the ratio of the number of intermediate- mass 
stars with m* = 4M - 5M for YS86 IMFs to that for Salpeter IMF is within a factor of 
3. If stars with mass of 5M are only formed, the number is larger by a factor of 15 than 
that for the Salpeter IMF with mi = 1M (see section 3). 



In this section, we show the evolution of bolometric luminosity, ionizing photon pro- 
duction rate, and emission line ratios of Population III star clusters with the YS86 IMFs. 
We adopt the 10 7 yr burst model as the star formation history. We calculate the spectral 
energy distributions (SEDs) of star clusters by using the evolutionary population synthesis 
code PEGASE2 (Fioc & Rocca-Volmerange 1997, 2000). Although the code can be used to 
calculate the photometric evolution of galaxies being the effect of the chemical evolution into 
account, we calculate that of a single metallicity population with Z = 0.0001; note that this 
metallicity is the poorest metallicity available in PEGASE2. The duration of the starburst 
is longer than the lifetime of the most massive star in our model (6OM ), 5 x 10 6 yr (Girardi 
et al. 1996). Since the lifetime is a half of the starburst duration, the metallicity of a half of 
stars formed in the starburst is Z < 10~ 10 . 



2.3. Galaxy Evolution Models 
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Figure 2 shows the SED evolution for each model: Salpeter IMF (solid lines), YS86a 
IMF (dotted lines) and YS86b IMF (dashed lines). We show the SEDs at 0.01 Gyr, 0.1 Gyr 
and 0.5 Gyr. Note that the ordinate is the flux in units of erg s™ 1 Mq 1 , corresponding to 
the case in which the total gas mass is 1M Q . The SEDs are mainly determined by the upper 
main sequence stars and are not sensitive to the IMF. During the starburst, the ionizing 
photons are radiated from massive main sequence stars whose effective temperature is lower 
than about 60,000 K (Girardi et al. 1996). Therefore there is a large break at 228 A; the 
photons with wavelength shorter than this are capable of ionizing the singly ionized helium. 
At an age of > 10 8 yr, hottest stars in the star cluster are not main sequence stars but post- 
asymptotic giant branch (AGB) stars whose progenitors are stars with a degenerate carbon 
core, e.g., intermediate- mass stars. The effective temperature of such post- AGB stars (i.e., 
PNNs) is about 200,000 K (Blocker 1995). Therefore the break of 228 A of such stars 
becomes to be small and thus the SED shape of far-UV is very hard. 

The global shape of SED depends slightly on the IMF (Figure 3). During a period of 
the starburst, between 4 x 10 6 yr and ~ 1 x 10 7 yr, colors (U — B, B — V, and V — K) 
become redder because of red supergiants. In this phase, the star cluster with the Salpeter 
IMF is reddest and that with the YS86a is bluest. This is because the number ratio of red 
supergiants to upper main sequence stars is highest for the Salpeter IMF (m -135 ) and lowest 
for the YS86a IMF (m~ 3 7 ). On the other hand, for ages longer than 10 8 yr, the star cluster 
with the Salpeter IMF is bluest and that with the YS86b is reddest. The reason for this is 
that the number ratio of redgiants to upper main sequence stars is lowest for the Salpeter 
IMF (m" L35 ) and highest for the YS86b IMF (m 3 - 5 ). 

Figure 4 shows the evolution of bolometric luminosity of a star cluster of 1M for 
various IMFs. The bolometric luminosity reaches its maximum at the end of starburst (10 7 
yr). The star cluster with the YS86b IMF is brightest while that with the YS86a IMF is 
faintest among the models. After that, following the death of massive stars, the bolometric 
luminosity of stellar systems decreases with time. The decreasing rate of the bolometric 
luminosity also depends on the functional form of IMF. Since the YS86a IMF has a peak at 
m m 4M whose lifetime is about 0.15 Gyr, the decreasing rate of the bolometric luminosity 
is small till the time, and becomes large after that. Since the lifetime of stars with 1M 
is about 5.6 Gyr, all the stars formed during the starburst have died till the time. As we 
show below, the dominant sources of the ionizing photons become the PNN at the age of ~ 
0.1 Gyr. The bolometric luminosity to mass ratios decrease from 66L /M to 3.2L /M 
during this period. 

In the upper panel of Figure 5, we show the evolution of production rate of photons 
capable of ionizing neutral hydrogen, A(H°), which is calculated by integrating the spectra 
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below 912 A, and that of photons capable of ionizing ionized helium, iV(He + ), which is 
calculated by integrating the spectra below 228 A. First, we mention how the production 
rate of the hydrogen ionizing photons evolves. As the bolometric luminosity increases, the 
ionizing photon production rate reaches its maximum at the end of starburst, because the 
lifetime of early B stars of ~ 1OM is longer than 10 7 yr. The evolution of iV(H°) depends 
on IMF. During the starburst, the ionizing photon production rate per unit mass is largest 
for the star cluster with the YS86b IMF but smallest for that with the YS86a IMF. In our 
calculations, we adopt the mass of most massive stars M u = 6OM whose lifetime is about 

5 x 10 6 yr. When the age of star clusters reaches the lifetime of the most massive stars, 
the number of hydrogen ionizing photons begins to decrease with time, following the death 
of stars radiating hydrogen ionizing photons. At the age of ~ 10 8 yr, the intermediate- 
mass stars evolve to PNNs and they become dominant ionization sources for hydrogen. In 
the star cluster with the Salpeter IMF, the number of PNNs hardly changes and thus the 
production rate of ionizing photons is kept nearly constant. On the other hand, for the star 
cluster whose IMF is either YS86a or YS86b, the number of PNNs and the production rate 
of ionizing photons decrease with decreasing number of low mass stars. After the age of ~ 
1 Gyr, the ionizing photon production rate is highest for the star cluster with the Salpeter 
IMF and smallest for that with the YS86b IMF since it is proportional to the number of 
PNNs. Next, we comment on the evolution of the production rate of ionized helium ionizing 
photons. The lifetime of stars with m > 10M & is shorter than a few xlO 7 yr. After the 
massive stars have died, the remaining main sequence stars cannot ionize the neutral helium. 
When the age of the star cluster become 1 x 10 8 yr, PNNs come into existence. The effective 
temperature of PNNs is very high enough to radiate He + ionizing photons. 

The lower panel of Figure 5 shows the evolution of the ratio of iV(He + ) to iV(H ) which 
indicates the hardness of UV spectra. During the starburst, the ratio is lower than 10~ 3 . 
However, after PNNs become the dominant source of UV radiation, the ratio exceeds 10~ 3 . 
The stellar mass of a progenitor of PNN decreases with time. Therefore, the effective tem- 
perature of PNN decreases and thus the ratio of iV(He + ) to iV(H°) also decreases (Vassiliadis 

6 Wood 1994). 

We investigate ultraviolet emission-line properties of the gaseous nebula photoionized 
by PNN clusters during the course of galaxy evolution. We use the photoionization code 
CLOUDY94 (Ferland 1997), which solves the equations of statistical and thermal equilibrium 
and produces a self-consistent model of the run of temperature as a function of depth into 
the nebula. Here we assume that a uniform density, dust-free gas cloud with plane-parallel 
geometry is ionized by various continuum sources. The parameters for the calculations are (1) 
the hydrogen density of the cloud, (2) the ionization parameter, U = (5(H°)(47rr 2 nHc) _1 , 
where Q(H°) is the number of ionizing photons, r is the distance from the ionizing source, 
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and c is the speed of light (Osterbrock 1989), (3) the SED of the ionizing radiation, and (4) 
the chemical composition. The SED of the ionizing radiation is derived as noted before using 
PEGASE2 (see Figure 2). We perform our photoionization calculations as a function of the 
ionization parameter between \ogU = —4 and —1 with a logarithmic interval of 0.5. We 
adopt nn = 10 2 cm -3 . We also adopt a metallicity of Z = 0.001 = 0.05Z Q . The calculations 
were stopped when the temperature fell to 3000 K, below which little optical emission is 
expected. 

For reference, we also calculate AGN photoionization models and shock-heating models 
in the following ways, a) AGN photoionization models are calculated by using CLOUDY94. 
We adopt that the input SED has a power-law form with a = —1.5 (/„ oc u a ). b) Shock- 
heating models are taken from results of Dopita & Sutherland (1996). Their models assume 
solar metallicities, shock velocities from 150 km s _1 to 500 km s" 1 , and the magnetic param- 
eter < Bjyjn < 4/iG cm 3 / 2 . The PNN photoionization model results are shown together 
with these models in Figure 6; the UV emission-line diagnostic diagrams proposed by Dc 
Breuck et al. (2000). 

Since high-energy photons supplied from the PNN clusters can cause high-ionization 
lines such as C IV, the UV emission-line properties of the PNN photoionization models are 
quite similar to those of the AGN photoionization. In particular, the model locus predicted 
for 10 7 yr-burst stellar populations with age of 0.5 Gyr cannot be distinguished from that 
for the AGN models with a power-law continuum with a = —1.5. This similarity has been 
already utilized to show that some LINERs (low-ionization nuclear emission-line regions: 
Heckman 1980) can be explained by photoionization by PNN clusters (Taniguchi, Shioya, 
& Murayama 2000). On the other hand, the shock- heating models appear to show different 
properties from those of the PNN and AGN models. 

Figure 7 shows the Lneii/^Boh Lya/He II, and N v/He II ratios as a function of logU 
for PNN models with age of 0.1 Gyr (open circles) and 0.5 Gyr (filled circles). These ratios 
hardly depend on the IMF but depend on the age. 

Recently, the evolutionary tracks of zero-metal stars have been published (e.g., Marigo 
et al. 2001). Based on these new models, Schaerer (2002) demonstrated that massive 
population III stars can radiate the He + ionizing photons because of both the high effective 
temperature (T cS > 80, 000 K) of zero-metal stars and the effect of the non-LTE stellar 
atmosphere (see also Tumlinson & Shull 2000). The small N(Re + )/N(B°) ratio during the 
starburst shown in Fig. 5 is attributed to the use of the population synthesis model with 
metallicity of Z = 0.0001 instead that with Z = 0. We note, however, H II regions around 
massive population III stars cannot radiate metal lines in which we are interested since the 
gas around massive Population III stars is still primordial. 
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2.4. Comparison with high-z star forming galaxies 

If the star formation in some forming galaxies takes place with YS86 IMF, such galaxies 
could experience a PNN-dominated phase ~ 10 8 yr after the initial starburst. Therefore it 
seems important to look for such galaxies at high redshift. As is shown in section 2.3, the SED 
of star clusters in PNN-dominated phase is so hard that the emission-line properties of the 
gaseous nebula photoionized by PNN clusters are indistinguishable from those photoionized 
by AGN. Considering the stellar evolution of intermediate mass stars, the PNN phase follows 
the post-AGB phase when stars eject a dusty gaseous envelope. It is therefore expected that 
the PNN cluster is surrounded by dusty gas (see the next section). Until the PNN-dominated 
phase begin, dust formed in the ejecta of Type II supernovae is also present in the interstellar 
medium (Todini & Ferrara 2001). Therefore, a part of ionizing photons could be absorbed 
by such dust grains. However, we neglect this effect in the following discussion for simplicity. 

The above expected features of the PNN star cluster are very interesting because the 
PNN star cluster may be observed as a dusty, ultraluminous infrared galaxy with high- 
ionization emission lines at high redshift. It is known that some high-z forming galaxy 
candidates share nearly the same properties; e.g., IRAS F 10214+4724 at z — 2.3 (Rowan- 
Robinson et al. 1991) and SMM 02399-0136 at z = 2.8 (Ivison et al. 1998). Their observed 
huge luminosities (L ho i > 1O 12 L ) have been often considered to be powered by massive 
stars. However, high-ionization emission lines such as C IV A1549 A and N v A1240 A are 
also observed in them. Since these emission lines cannot be produced by photoionization by 
ordinary massive stars, it has been often considered that active galactic nuclei (AGNs) also 
reside in them (Ivison et al. 1998; Serjeant et al. 1998). In particular, the presence of hidden 
broad emission lines reinforces this interpretation for IRAS F10214+4724 (Goodrich et al. 
1996). Spectropolarimetry of SMM J02399— 0136 reveals that the polarization behavior is 
similar to BAL quasars (Vernet & Cimatti 2001). Even in the local universe, bursts of star 
formation (starbursts) are sometimes associated with galaxies hosting AGNs and thus the 
simultaneous presence of both starburst and AGNs may not be unusual (e.g., Taniguchi 
1999; Mouri & Taniguchi 2001; Storchi-Bergmann et al. 2001). However, the evidence for 
AGNs in them has been obtained often by rest-frame ultraviolet and optical emission lines 
which are subject to misunderstanding in some cases (Taniguchi et al. 1999; Lutz, Veilluex, 
& Genzel 1999). Since the effective temperatures of PNNs exceed 10 5 K, high energy photons 
are much more numerous than those from ordinary massive stars. Here we examine whether 
or not the PNN cluster is responsible for the formation of high-ionization emission lines. 

In Figure 6, we plot the observed positions of IRAS F10214+4714 (Serjeant et al. 1998) 
and the high-z powerful radio galaxies studied by De Breuck et al. (2000). Although the 
high-z powerful radio galaxies have an AGN, at least some of them are considered to be very 
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young galaxies (e.g., Eales et al. 1993). We also plot available data of planetary nebulae 
(Feibelman 2000). It is found that the emission line ratios of IRAS F10214+4724 as well 
as the high-z powerful radio galaxies in these UV diagnostics can be explained either by 
the AGN models with —2 < loglf < —1.5 or by the PNN model with an age of 0.5 Gyr. 
The similar behaviors between the PNN and AGN photoionization models unable us to 
distinguish the two models unfortunately. 

It is noted that both the AGN and PNN models appear to underpredict the C n]/C ill 
ratio as shown in the right-lower panel in Figure 6. This discrepancy suggests that the ob- 
served UV emission-line ratios of IRAS F10214+4714 and the high- ,2 powerful radio galaxies 
can be explained either by a combination between the shock heating and the AGN photoion- 
ization or by a combination between the shock heating and the PNN photoionization. 

Next, we study whether or not the star clusters dominated by PNNs can explain the 
luminosity of high-z forming galaxy candidates. As shown in section 2.2 the bolometric 
luminosity to mass ratio changes from 66L /M© to 3.2L /M . Therefore, the mass of stellar 
system needed to reproduce the bolometric luminosity of IRAS F10214+4724 (~ 1O 13 L ), 
is 10 11 — 1O 12 M . This value is comparable to the stellar mass of massive galaxies. It is 
impossible to radiate the huge luminosity from only PNNs of first generation stars for this 
galaxy. There are two reasons why it is hard to explain the huge luminosity by a lot of 
PNNs. As shown in section 2, even for the YS86's IMF the number of intermediate-mass 
stars (e.g., 4M < M* < 5M ) is only a few time larger than that for the Salpeter IMF. 
The other reason is related to the difference of the lifetime of progenitor stars. The lifetime 
of stars with mass of 4M is 1.54 x 10 8 yr and that of 5M is 9.9 x 10 7 yr (Girardi et al. 
1996). Therefore, the mass range of stars which evolve to PNNs at the same time is very 
small and as a result the number of those stars is very small. To solve this discrepancy, we 
consider an extreme model in the next section. 



3. AN EXTREME MODEL 
3.1. Scenario 

Here, we consider an extreme case that the formation of intermediate-mass stars occurs 
during a period of tsf- Most luminous forming galaxy candidates at high redshift contain 
a lot of molecular gas (e.g., M gas ~ 1O U M ) and their inferred star formation rates (SFR) 
are of the order of ~ 1O 3 M yr -1 (Frayer et al. 1998, 1999). Therefore, we assume that a 
galaxy has a gas mass of 1O 11 M and the initial star formation occurs at a rate of SFR = 
Msf = 1O 3 M yr -1 . We adopt a duration of the star formation tsf = 10 7 yr. The gas mass 
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converted to stars is estimated as M S f = Msf x tsf = 1O 1O M . Thus in this model, the star 
formation efficiency is t]sf = 10%. The star formation may take place in a central region 
(say, within a radius of ~ 1 kpc) of galaxies. We assume for simplicity that all stars formed 
in this initial star formation have a mass of 5 M (i.e., a typical intermediate-mass star). 
Therefore, the SFR adopted above means that 200 stars with m* = 5M are formed in a 
year; AT* = Msf/^* = 200 stars yr -1 . 

Here we consider the evolution of a cluster of such stars. The bolometric luminosity of 
each star in the main sequence phase is L* « 1 x 1O 3 L and the age of the star is r age ~ 1 x 10 8 
years. The total number of stars formed in the duration is iV* (total) = iV* x tsf = 2 x 10 9 
stars. Therefore, the total luminosity of the star cluster increases during the star formation 
period (10 7 years) at a rate of iV* x 2 x 1O 5 L yr -1 and reach the maximum bolometric 
luminosity Lboi(max) ps 2 x 1O 12 L at t — 10 7 years. This maximum luminosity keeps as 
long as the stars are alive; i.e., up to t = r agc = 10 8 years. Shortly after this age, the stars 
evolve from the main sequence to PNNs through AGB. 

Since one PNN has a bolometric luminosity of L PNN ks 1O 4 L (e.g., Vassiliadis & Wood 
1994), the total PNN luminosity amounts to L PNN (total) = iV* r PNN L PNN = 2 x 1O 11 L 
where t P nn ~ 10 5 yr is the lifetime of the PNN phase. Nevertheless this short lifetime, note 
that the PNN phase in the star cluster continues during nearly the same period as that of 
star formation; i.e., 10 7 yr. In the PNN phase, the total luminosity of the star cluster is 
mostly dominated by unevolved main sequence stars while the ionizing photons are supplied 
mainly by PNNs. 

Another important property of the PNN phase is that a lot of gas is ejected from the 
stars. The core mass of each star is estimated to be m c ~ O.6M . Therefore, each star ejects 
a metal-enriched gas of m gaSjP N = — m c ~ 5M — O.6M ~ 4.4M . The total ejected gas 
mass amounts to M gas P N = m gaS) PNiV* (total) ~ 8.8 x 1O 9 M . Since the gas-to-dust mass 
ratio in such gas may be ~ 100, the star cluster can produce a lot of dust grains with a mass 
of M dust ~ 1O 8 M . 

In summary, the PNN phase of the star cluster has the following properties; 1) the 
bolometric luminosity is L ho \ ~ 1O 12 L ; 2) the ionizing luminosity is ~ 2 x 1O 11 L ; 3) the 
ionizing continuum is so hard that it creates highly-ionized regions around the star cluster; 
4) the dust mass is M dust ~ 1O 8 M ; 5) the heavy element abundance is increased up to 
Z ~ M dust /M gas ~ 0.001, and 6) this phase lasts for xs 10 7 yr given that the duration of star 
formation is 10 7 yr. 
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3.2. Emission-line properties of the PNN cluster 

As shown in section 2.3, the PNN cluster is capable of the formation of the high- 
ionization lines such as C IV A1549 and C ill] A1909 emission lines. For example, the 
observed C iv/He II ratios are 5.7 and 1.9 for SMM 02399-0136 and IRAS F10214+4724, 
respectively (Ivison et al. 1998; Serjeant et al. 1998). However, our models cannot explain 
the observed higher ratios of N vA1240/He II; 8.5 for SMM 02399-0136 (Ivison et al. 1998) 
and 1.4 for IRAS F10214+4724 (Serjeant et al. 1998) (see Figure 7). In our model, we 
have assumed that each heavy element abundance is five percent of its solar value. If only 
the helium burning occurs in the core of PNN, carbon could be enriched. However, if the 
hydrogen burning also occurs, nitrogen could be enriched because nitrogen is formed as the 
end product of CNO cycle in the core. If such nitrogen enrichment occurs, the observed 
strong N v emission may be explained. It is worthwhile noting that the overabundance 
of nitrogen is also suggested for the broad emission-line gas in high- ,2 quasars (Hamann & 
Ferland 1993, 1999). 

3.3. Further Evolution of the Star Cluster 

We consider the final fate of the star cluster in our model given in section 3.1. The 
intermediate-mass stars formed in the initial starburst will evolve to white dwarf (WD) 
stars with a mass of m WD = m c m O.6M through the phases of AGB and planetary 
nebula. Since the evolved WD stars have a luminosity of ~ 1O~ 4 L at their ages of 10 10 
yr, the star cluster formed in the starburst becomes finally very faint; i.e, Lwd (total) = 
Lwd-'V* (total) ~ 2 x 1O 5 L . Since there is no such candidate in the local universe, it is 
expected that further star formation could have occurred in galaxies even if starbursts with 
only intermediate-mass stars occurred in the galaxies. After the death of PNN stars, dust 
grains with a mass of ~ 1O 8 M are left in a galaxy. Since the galaxy has the gas with a 
mass of 1O 11 M , the heavy element abundance increases at a level of Z ~ 0.001 after the 
PNN phase. Therefore, so-called Population II stars would form in the gas clouds polluted 
chemically by the intermediate-mass stars. In other word, galaxies with a PNN cluster may 
evolve to ordinary galaxies observed in the present-day universe. 

3.4. Final Remarks 

What we would like to propose in this section is not that some particular high- ,2 dusty 
galaxies are indeed the PNN galaxies but that candidates of the extreme PNN galaxies share 
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nearly the same properties of the observed high- z dusty galaxies. Recent deep surveys both 
at FIR and submillimeter wavelengths have discovered a significant number of hyperluminous 
dusty galaxies at high redshift (Rowan-Robinson et al. 1991; Kawara et al. 1998; Puget et 
al. 1999; Smail, Ivison, & Blain 1997; Hughes et al. 1998; Barger et al. 1998; Eales et al. 
1999; Rowan- Robinson 2000). If some proto-galactic systems tend to form intermediate-mass 
stars preferentially as Population III stars, it is possible that candidates of the PNN galaxies 
may be found in some high-z dusty galaxies. If present, such galaxies will be detected in 
future sensitive search for dusty galaxies made by Atacama Large Millimeter Array (ALMA; 
Kurz & Shaver 1999). Deep optical spectroscopy using new 30-m class telescopes will be 
also necessary to investigate detailed rest-frame UV spectra of such high-z dusty galaxies. 
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Table 1. The number of intermediate-mass stars with 4M < < 5M . 



IMF 


iV(4M < ///. < 5M ) 


Salpeter (M, = O.1M ) 


0.0052 


Salpeter (M, = 1M ) 


0.0136 


YS86 (77 = 3)(YS86a) 


0.0323 


YS86 (77 = 1.5)(YS86b) 


0.0149 


5-function (5M ) 


0.2 
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Fig. 1. — Comparisons of YS86 IMFs (dotted line: YS86a, dashed line: YS86b) with the 
Salpeter IMF (solid line). 
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Fig. 2. — The spectral energy distributions of star clusters with age of 0.01 Gyr (upper lines), 
0.1 Gyr (middle lines) and 0.5 Gyr (lower lines) for Salpeter IMF (solid lines), YS86a IMF 
(dotted lines) and YS86b IMF (dashed lines). 
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Fig. 3. — Evolution of U — B, B — V and V — K colors of 10 7 yr starburst model with 
Z=0.0001 for various IMFs: Salpeter (solid lines), YS86a (dotted lines) and YS86b (dashed 
lines). 
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Fig. 4. — Evolution of bolometric luminosity of star clusters of mass of 1M for Salpeter 
IMF (solid line), YS86a IMF (dotted line) and YS86b IMF (dashed line). 
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Age (yr) 

Fig. 5. — (upper) Evolution of ionizing photon production rate with wavelength below 912 
A(7V(H°), thick lines) and that with wavelength below 228 A(iV(He + ), thin lines) for var- 
ious IMFs: Salpeter (solid lines), YS86a (dotted lines) and YS86b (dashed lines), (lower) 
Evolution of the ratio of iV(He + ) to iV(H°) for various IMFs: Salpeter (solid lines), YS86a 
(dotted lines) and YS86b (dashed lines). 
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Fig. 6. — Ultraviolet emission line diagnostics adopted in De Breuck et al. (2000) to discuss 
an excitation mechanism. Black lines with black marks are PNN models at 0.1 Gyr (open 
circles) and 0.5 Gyr (filled circles) for various IMF: Salpeter (solid lines), YS86a (dotted 
lines) and YS86b (dashed lines). Blue lines with blue marks are AGN models, power law 
with a = —1.5 (f u oc u a ) for various density: n e = 10 2 cm -3 (solid lines), n e = 10 3 cm -3 
(dotted lines), and n e = 10 4 cm -3 (dashed lines). Red lines with red marks are shock model 
calculated by Dopita Sz Sutherland (1996) for various magnetic parameter: B / ' yfn = (solid 
line), 1 (dotted line), 2 (short dashed line) and 4 (long dashed line) \i G cm 3 / 2 . Green marks 
show observational points: IRAS F10214+4724 (circle), high-z radio galaxies (crosses, De 
Breuk et al. 2000), and planetary nebulae (filled triangle, Feibelman 2000). The dereddening 
vector (Ay = 1) is shown in each panel. 
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Fig. 7. — Dependence of log(LH e ii/ -^Boi) , log Lya /Hell, and log Nv/Hell on log U in our PNN 
models at 0.1 Gyr (open circles) and 0.5 Gyr (filled circles) for various IMF: Salpeter (solid 
lines), YS86a (dotted lines) and YS86b (dashed lines) 



